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We study the gravitational collapse in (n+2)-D quasi-spherical Szekeres space-time (which possess 
no killing vectors) with dust as the matter distribution. Instead of choosing the radial coordinate 'r' 
as the initial value for the scale factor R, we consider a power function of r as the initial scale for the 
radius R. We examine the influence of initial data on the formation of singularity in gravitational 
collapse. 

PACS numbers: 0420D, 0420J, 0470B 

Over the last two decades gravitational collapse in spherical space-time (TBL model) has 
been studied extensively with dust as the matter content [1-10]. A general conclusion from 
these studies is that a central curvature singularity may be naked but its local or global 
visibility depends on the choice of initial data. If the regular initial density profile falls 
off rapidly having a maximum value at the centre then it is not possible to have naked 
singularity above five dimensional space-time [8-10]. 

Also from the recent past, attention has been given to study non-spherical collapse 
[11-20]. Most of these studies deal collapse numerically [12-15] with a few analytical 
works [16] (for quasi-spherical gravitational collapse, see ref. [17-20]). These are mainly 
concerned with special shape of the gravitating body. The present work examines the 
role of initial data in the formation of gravitational collapse in (n + 2)-D Szekeres space-time. 

The metric ansatz for {n + 2)D Szekeres space-time is [21, 22] 



ds^ = df - e"^dr^ - 6^/3 ^ dx^^ (1) 

i=l 

where a and (3 are functions of all the (n -t- 2) space-time co-ordinates. But if we assume 
that /3' ^ then for inhomogeneous dust model the solution is [22] 

e'' = i?(t,r) e^('^'^i (2) 



e" = (3) 



i=l i=l 



n F' + {n + l)Fiy' 



p{t,r,xi,...,Xn)- ^ i?»(i?' + i?,.') 



(6) 
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-AAC = -1 



(7) 



i=l 



where A{r), Bi{r) and C(r) are arbitrary functions of r satisfying (7) and also in the 
expression (5), /(r) and F{r) are arbitrary functions of r alone. 

A shell focusing singularity on a shell of dust occurs when it collapses at or starts ex- 
panding from the centre of matter distribution. Wc shall consider only the central shell 
focusing singularity (i.e., i? = or /3 = — oo) for marginally bound case only (i.e., /(r) = 0). 
Suppose t = ti he the initial hypersurface from which the collapse develops. For initial 
data we assume that R{ti,r) is a monotone increasing function of r. So without any loss of 
generality, it is possible to make an arbitrary relabeling of the dust shells by r ^ g{r) such 
that we can choose 



R{ti,r) = Ror'^, {Rq > 0, k are constants) 
Hence solving equation (5) using the initial condition (8) we get 

i+i (n+i)fc n+1 



R 



Rci 



r 2 



(8) 



(9) 



The regularity of the metric co-efficients on the initial hypersurface restricts A: > 1 . Further 
for the regularity of the initial density profile 



. , s_ri F' + {n + l)Fv' 

PiV ' "^1' "^nj — PV'ii'^i-^li •••I'^n) — „ r>n-\ 



2 i?J+V("+i)'=-i(fc + r!/') 
we can write the following series expansion for Pi(r), F{r) and v'{r) near r = as [22] 



(10) 



3=0 



(11) 



F{r) = r("+i)'=+-'' 



(12) 



and 



u'{r) = J2 ^0 r\ (I'-i + A; > 0) 



(13) 



Also using equations (11)- (13) in (10) we have the following relations among the different 
co-efBcients 



2 P^=2l-+1+A + . 
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-(n+l) 
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n+ 1 + 
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R, 



-(n+l) 
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-(n+l) 
•0 ' 

(15) 



according as > —k or i/ .^ = —k. 

The singularity curve t = ts (r) for the shell focusing singularity is characterized by 



i?(t,(r),r) =0 

So the time for central shell focusing singularity is given by 



(16) 



2R ^ 
(n + l)V-Pb 

Further, if tahir) is the instant for the formation of apparent horizon then we have 



(17) 



which gives 



R''-\tahir),r) = F{r) 



(18) 



tah{r) -to = - 



Rn 



(n + l)Fo'/' 



n + l 



{n + l)k 1 

r + 



n-lFo 



r "-1 



+1 



+ . 



(19) 



The above expression shows the time difference between the formation of trapped surface 
at a distance r and the time of singularity formation at r = (central singularity). Hence 
the necessary condition that an observer at a distance r will observe the central singularity 
(at least locally) is tahi''') > to (for details see Ref. [6—9, 24]). 

From physical consideration it is reasonable to assume that the initial density pi{r) is 
maximum at the centre r — 0. This implies that the first non- vanishing term after po in the 
scries expansion (see eq. (11)) for Pi{r) should be negative. Further, one may assume that 
Pi{r) should vanish at r = but is negative in the neighbouring region or more generally, it 
may be assumed that Pi{r) = p'l{r) = ... = p'l^^ir) — and pf{r) < (s > 2). Now using 
the relations (14) or (15) among the co-efBcients we have 
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I. when > —k: 

Here pi < implies Fi < 0. Also in general, pi = P2 = ■■■ = Ps-i = and Ps < implies 
Fi=F2 = ... = Fs-i = and < with s > 2. 

II. when J/ J = —k: 

Here pi < does not imply Fi < 0. Also more generally, pi = p2 = ... = Ps-i = and 
Ps < may have one possible solution as Fi = F2 = ... = = and Fg+i < with s > 2. 

In particular, if we assume the initial density to have a maximum value at r = and 
falls off rapidly near r = then we have pi = 0, p2 < near r = 0. Then in case I, we have 
Fi = 0, F2 < while in case II, we may take Fi = F2 = 0, F3 < near r = 0. 

Therefore in the present problem we have the following possibilities for naked singularity: 

(a) If Fi < then naked singularity may appear in any dimension (> 4) for v_-^ > —k 
and fc > 1. 

{h) In general if we choose F\ = 0, F2 = 0, ...,Fj_i = and F, < (i > 2) then for 
formation of naked singularity, 'n' is restricted by the inequality 



(20) 

with max{l, |) < /c < i, > —k, k > 1. 

Here [x] stands for the greatest integer in x. From the inequality (20) wo note that if 
'fc' is very close to H' (but less than i ) then 'n' can take larger values than 2 i.e., naked 
singularity may appear in much larger dimension compare to the usual four dimension. The 
following table shows some possible values of n for different values of i and k (with k < i) 
from the inequality (20). 



2 < n < 



i + k 



TABLE: Some possible values of n for different values of i and k. 



i — > 
ki 


2 


3 


4 


5 


3 
2 


n = 2,3,4,5,6,7 


n = 2,3 


n = 2 




21 
11 


2 < n < 41 


n = 2,3,4 


n = 2 




35 
12 




2 < n < 71 


n = 2,3,4,5,6 


n = 2,3 



For > —k we have from equation (14) Fi = 0, F2 < and hence for naked singularity 



2 < n < 



'2 + k' 
2-k 



with 1 < fc < 2, while for 
restricted by the inequality 



-k we choose Fi = 0, F2 = 0, F3 < and so n will be 



5 



2<n< 



"3 + fc" 
3- fc 



with 1 < A; < 3. 



We note that for k — 1 the possible values of n are 2 and 3 for > — 1 and n = 2 for 
= — 1 i.e., naked singularity is possible only for four and five dimensions which we have 
shown in earlier works (see ref. [8-10, 20]). 

Now we shall examine the nature of singularity by studying the outgoing radial null 
geodesic (ORNG) originated from the central shell focusing singularity. Let us start with 
the assumption that it is possible to have one or more such geodesies and we choose the 
form of ORNG in power series as [23, 24] 



t = to+ar^, (21) 

upto leading order near r = in t-r plane with a > 0, ^ > as constants. Using equation 
(16) and (17) the singularity curve can be written as (near r = 0) 

n+1 

tsir)^to- ^"'^\,^ r^ (22) 

where m > 1 is an integer and is the first non- vanishing term beyond _Fo . As for naked 
singularity we have t < ts{r) so comparing (22) with (21) for ORNG the restrictions on the 
two parameters ^ and a as 

n+1 

e>m and a< ^"'^o ' , (23) 

(n + l)Fo'/' 

Moreover, from the metric (1) we can write for ORNG 



^=R' + Ru' (24) 

ar 

We shall now examine the feasibility of the null geodesic starting from the singularity 
with the above restrictions for the following two cases namely, ^ > m and ^ = m. 

When ^ > TO then near r = the solution for R in (9) (choosing ti = 0) simplifies to 



R = Ro{-^] r^^^' (25) 



2Fo 

Now combining (21) and (25) in equation (24) we get (upto leading order in r) 



««/- = ^(._. + . + ^)(-^) (2e) 



which implies 
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Now if < m then n and k are bounded by the same inequalities as in (20) (with i = m) 

for the formation of naked singularity while for k > m, there will be no upper bound of n. 
Furthermore, we note that from equation (27), ^ > and a > as ^ > —k. Thus we 
have the same conclusion as before (in case (b)) and it is possible to have consistent ORNG 
originated from the singularity. 

On the other hand for ^ = m, we have from equation (24) using the solutions for R and 
v as before, it is possible to have naked singularity if 



m + k , 

n= -, max(l,—j<k<m (28) 

m — k 3 



and 



a = — - 



1 / n+ 1 



m 



P n+l \ 1+ 

-Kn 



2 2Fn ° 



(29)" 



Now we shall examine whether the restriction (23) for 'a' is consistent with the expression 
'a' in equation (29). In fact equation (29) takes the form 



2^6m = -[-(n+l)6-C]^ [{u + m)C + [n + l){v + k)b] (30) 
with the transformation 

a = bFo-\ Fm= ^ . (31) 

Since equation (30) is a real valued equation of 6, so we must have 

(n+l)6 + C<0 , 

which using (31) gives us the restriction on a in equation (23). Hence the geodesic (21) will 
have consistent solution for 'a' and So the above conclusion regarding the formation of 
naked singularity is justified. Further, introducing the variable 'i?^' by the relation 

0=-(n+l)6-C (32) 

we have seen from equation (30) 

4A;"+V""'(C + (PT^^ = [2fcC - (n - + fc)</>]"+^ (33) 

with the restriction < ^ < — C and 'a' will satisfy the inequality in equation (23). 

Further we observe that for k > m, n can not have any positive integral solution and 

hence naked singularity is not possible. Also this case (i.e., ^ = m) has no analogue with 
our previous result by comparing tahiT) with to- Lastly for any fixed m, k can only take 
those values within the limits in equation (28) which makes n, a positive integer (> 2). It is 

to be noted that if we consider next order term in the geodesic eq. (21) then it is possible 
to have a family of ORNG originated from the central singularity (see Ref. [24] ) . 

Thus the formation of naked singularity strongly depends on the nature of the initial 
density and also on the choice of the parameter k involved in the initial choice of the scale 
factor R. Therefore we conclude that formation of naked singularity strongly depends on 
the choice of initial data for the physical parameters as well as for geometric quantity. For 
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future work, it will be interesting to study the dominance of initial data for physical pa- 
rameters over that for geometric quantities and vice versa for formation of naked singularity. 

Acknowledgement : 

One of the authors (U.D) thanks CSIR (Govt, of India) for the award of a Senior 
Research Fellowship. 

References: 

[I] A. Ori and T. Piran, Phys. Rev. Lett. 59 2137 (1987). 
[2] F. J. Tipler, Phys. Lett. A 64 8 (1987). 

[3] F. C. Mcna, R. Tavakol and R S. Joshi, Phys. Rev. D 62 044001 (2000). 

[4] P.S. Joshi, Global Aspects in Gravitation and Gosmology (Oxford Univ. Press, Oxford, 

1993). 

[5] P. S. Joshi and I. H. Dwivedi, Commun. Math. Phys. 166 117 (1994); Class. Quantum 
Grav. 16 41 (1999). 

[6] D. Christodoulou, Commun. Math. Phys. 93 171 (1984); R. P. A. C. Newman, Class. 
Quantum. Grav. 3 527 (1986). 

[7] P.S. Joshi, N. Dadhich and R. Maartens, Phys. Rev. D 65 101501(i?)(2002). 

[8] A. Banerjee, U. Debnath and S. Chakraborty, 5r-gc/0211G99 (2002) (accepted in Int. J. 

Mod. Phys. D) . 

[9] U. Debnath and S. Chakraborty, gr-qc/ (2002) . 

[10] R. Goswami and P.S. Joshi, ,gr-gc/021 12097 (2002) . 

[II] C. Barrabes, W. Israel and P. S. LeteUer, Phys. Lett. A 160 41 (1991); M. A. Pelath, 
K. P. Tod and R. M. Wald, Class. Quantum Grav. 15 3917 (1998); F. Echeverria, Phys. 
Rev. D 47 2271 (1993); T. A. Apostolatos and K. S. Thorne, Phys. Rev. D 46 2435 (1992). 
[12] S. L. Shapiro and S. A. Teukolsky, Phys. Rev. Lett. 66 994 (1991); Phys. Rev. D 45 
2006 (1992). 

[13] T. Nakamura, M. Shibata and K.I. Nakao, Prog. Theor. Phys. 89 821 (1993) . 

[14] T. Harada, H. Iguchi and K.I. Nakao, Phys. Rev. D 58 041502 (1998) . 

[15] H. Iguchi, T. Harada and K.I. Nakao, Prog. Theor. Phys. 101 1235 (1999); Prog. 

Theor. Phys. 103 53 (2000). 

[16] K. S. Thorne, 1972 in Magic Without Magic: John Archibald Wheeler, Ed. Klauder J 

(San Francisco: W. H. Freeman and Co.). 
[17] P. Szekeres, Phys. Rev. D 12 2941 (1975). 

[18] P. S. Joshi and A. Krolak, Class. Quantum Grav. 13 3069 (1996); S. S. Deshingkar, S. 
Jhingan and P. S. Joshi, Gen. Rel. Grav. 30 1477 (1998). 

[19] S. M. C. V. Goncalves, Class. Quantum Grav. 18 4517 (2001); Phys. Rev. D 65 
084045 (2002). 

[20] U. Debnath, S. Chakraborty and J. D. Barrow, .9r-gc/0305075 (2003). 

[21] P. Szekeres, Commun. Math. Phys. 41 55 (1975). 

[22] S. Chakraborty and U. Debnath , ,9r-gc/0304072 (2003). 

[23] S. Barve, T. P. Singh, C. Vaz and L. Witten, Class. Quantum Grav. 16 1727 (1999). 
[24] U. Debnath and S. Chakraborty, math-ph/mQ7Q2A. 



